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1. Introduction 

Deep-inelastic scattering (DIS), measured in fixed-target experiments and at HERA, provides core 
constraints on the parton distributions for the LHC. For some crucial processes, such as gauge- 
boson and Higgs production, these distributions are required at the next-to-next-to-leading order 
(NNLO) of perturbative QCD. Consequently coefficient functions at this accuracy are needed also 
for the extraction of the parton densities from (mainly) the structure function F2{x,Q-^) in DIS. 

For the massless case, these quantities have been known for a long time |jT]]. However, a 
considerable part of F2 at small Bjorken-x is due to the production of charm quarks which is domi- 
nated by the photon-gluon fusion process Y*g ccX. The NLO coefficient functions for have 
been obtained in a semi-analytic manner iQ] ; the results are often used via the parametrizations of 
Ref. [^] (for minor corrections see Ref. [Q]). The corresponding NNLO corrections are not known. 
Fully analytic NLO results have obtained in the asymptotic limit m^/Q^ — )• [^. Recently these 
calculations have been extended to NNLO for the lowest even-integer Mellin moments [^]. 

It has been known for a long time, see, e.g., Refs. |^], that, at not too large values of Q^, the 
convolution of the coefficient function for Fj and the gluon density is dominated by rather low 
partonic of-mass energies (CM). Hence the NNLO predictions of the threshold resummation M, 



can provide useful information on the dominant contribution to Fj. Previously the first two [IC] 
and three [jll]] highest threshold logarithms have been determined at this and all higher orders. 

In this contribution we employ recent developments concerning the structure of massive- 
particle amplitudes and the description of heavy-quark production in hadronic collisions [|^, [I^ to 
extend those results to four logarithms, i.e., we are now able to derive all threshold-enhanced terms 



at NNLO. We also include a brief update of the results of Ref. [jlOp for the transverse momentum 



distributions, calculated at NLO in Refs. [14], using a modern set of parton distributions [15] 



2. Threshold resummation of the gluon coefficient function for F2 

The heavy-quark coefficient functions for F2 are usually expressed in terms of the variables 

= ^ , and 77 = 1 or jS = a/1 -P with p = , (2.1) 

m-^ p s 

where s is the CM energy, m the mass of the heavy quark, and jS the relative velocity of the heavy- 
quark pair. In terms of the threshold limit, p corresponds to the Bjorken variable x in massless 
DIS. Hence the dominant gluon coefficient function receives a double-logarithmic higher-order 
enhancement at j8 ^ 1 . The resummation of these logarithms is performed in terms of the Mellin 
variable N conjugate to p. Up to terms suppressed by powers of 1 /N, the coefficient function reads 

C2,AoCs,N) = 4^hN)-goias,N)-cxp[G{asMN)] . (2.2) 



Here is the lowest-order coefficient function (see, e.g., Ref. and gQ{(Xs,N) a matching 
coefficient. Its dependence on A'^, absent in the massless case, is due to Coulomb terms which aie 
enhanced by a factor l/p (see below). The resummation exponent G is of the standard form 



G= dz'— / 4-A,(a3(^')) +£>r?^cf(«s(W[l-z]2)) 

JO 1 - Z L-'M^ " 



(2.3) 
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The first term resums the coUinear gluons emitted off the initial gluon, the corresponding 'cusp 



anomalous dimension' A„ is known to order [16|. The second term collects soft and final-state 



emissions. Following the methods of Refs. [ |12| , [T^ ] we find 

Dyg^cc = ^/^Dgg^Higgs + Dqq , (2.4) 



where the latter heavy-quark coefficient is known to order [13] (obviously only the colour-octet 
result is required in the present case), and the former even to order [jl7|]. 

The above information is sufficient to predict the highest four powers of InN at all orders in as 
(cf., e.g., Ref. [[l8|]), provided that the matching function go is known at NLO. go is of the form 

goias,N) = gliia,) ■ gl^ia^N) . (2.5) 

The Coulomb contribution gQ can be determined by Mellin transforming the partonic cross section 



in non-relativistic QCD, calculated for the colour-singlet case to NNLO in Ref. [|19|]. The required 
octet results are obtained by the colour-factor replacement Cf — )• Cf — Ca/2. The NLO contribu- 
tion to the A/^-independent hard matching constant g^ had not been determined before this research. 
We have extracted this coefficient - which will be presented elsewhere [^] - analytically by inte- 
grating the intermediate results of Ref. [Q] (distributed as a FORTRAN program), and checked our 
result numerically, for some relevant values of , using the parametrization of Ref. 

3. Threshold approximation to the NNLO coefficient function 

The above A^-space results can be readily expanded in oCs and then Mellin inverted using, e.g., 
App. A of Ref. [ p| ] and the fact that the leading-order coefficient function is hnear in j8 near 
threshold (we normahze the coefficient functions as in Refs. [§,§]), 

cfli^,l5) = 7rr^j8(l + ^/4)-i + ^(j83) . (3.1) 



At NLO one thus recovers the threshold expansion (with Ty = 1/2, Ca = 3 and Cf = 4/3 in QCD) 

(0) 2 

= ^{4CAlnHsp^)-20CAln{SI5^) + co{^)+{2CF-CA)j 

+ ln4[-4CAln(4j82) + c-o(^)]+^(^2a ^32) 

The logarithmic and l/j3 contributions have first been given in Ref. [^]. The scale term Cq(i^) is 
fixed by renormalization-group constraints and reads 

c-o(^) = 4CA(2 + ln(l + ^/4)) -4/3r/ , (3.3) 

where the final term arises from the transformation of to the standard MS scheme [^] which was 
not performed in Ref. The corresponding scale-independent contribution Cq(i^) is not available 



in the literature yet, the full result will be presented in Ref. [20]. Here we can, for brevity, only 



provide its numerical values at the two scales used in our illustrations below, 

Co(1.956) = 88.28 , Co(19.56) = 70.23 . (3.4) 
For F2, hence nj: = 3 light flavours, our corresponding new NNLO results are numerically given by 
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JO) 



cg(^,i8) ^ ^{ln4^1152-W/3(l545. + 1152L) 



+ ln2/3 (^-3570.+48co(i^) + (118.0 + 48co(^))L + 288L2-167r2j8-' 
+ lnj8(^2403.-20.19co(^) + (2223.-20.19co(i§)-24co(^))L 

+ (291.3 -24co(^))l2 + 7r2^-' [2.910 + 8L]) + 0{p-^)^ (3.5) 

with L = ln(/i^/m^), where the coefficients with a decimal point are approximate. In addition to 
the terms given here, also the non-logarithmic 1 /j8 Coulomb contributions are now known. 
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Figure 1: Successive approximations to the NNLO gluon coefficient function for F2 in terms of threshold 
logarithms and 1 /j3 Coulomb contributions at two typical scales for a charm pole-mass m = 1 .43 GeV. 



The threshold expansion ( |3.5| ) of the NNLO coefficient function is shown in Fig. 1 for a stan- 
dard choice of the renormalization/factorization scale jJ.. Keeping only the highest two logarithms 
is obviously insufficient. The new InjS contribution is rather small at the lower, but definitely 
relevant at the higher scale, while the non-logarithmic NNLO Coulomb terms are small in both 
cases. The resulting estimates for the NNLO corrections to are illustrated in Fig. 2. In the 
region 10^^ ^ x < 10^^ these amount to no more than about 5 — 10% at = 40GeV'^, but reach 
15 — 30% at 2^ = 4 GeV^, which the largest effects occurring at the upper end of the above A;-range. 
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Figure 2: NLO and threshold-estimated NNLO results for the charm contribution to the structure function 
F2 using the respective parton distributions and strong coupling constants of Ref. [[l5[] with mc = 1 .43 GeV. 
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4. The pj^-differential charm structure function 



Experimentally the inclusive structure function is determined via (theory-dependent) extrapola- 
tions of more differential cross sections. As an example we consider the /j^-unintegrated structure 
function dF2/dpj, calculated at NLO in Refs. [14]. First NNLO estimates based on the next-to- 
leading logarithmic (NLL) threshold resummation were derived in Ref. [|lO[]. In Fig. 3 we present 



an update of these predictions, using an independent code and up-to-date parton densities [15] 



-1 

10 




/ dp^ (GeV ') 






x = 0.001, Q-= lOGeV^ 


10 " 




NNLO(NLL) 






NLO(fuIl) 






. . NLO(NLL) 


-3 

10 


' ^i^ = Q^ + 4m^ + 4p^ 






ABKM(09) PDFs 







dF^/dp^ (GeV ') 




x = 0.01,Q-= lOGeV- 




NNLO(NLL) 




NLO(full) 




\ NLO(NLL) 


= + 4 + 4 




ABKM(09) PDFs 






Al^ 



PT(GeV) 



4 5 
Pt (GeV) 



Figure 3: NNLO estimates for the /:ij,-unintegrated charm structure function F2 for two typical values of x. 
At NLO the results for the NLL expanded coefficient function are compared compared to the exact values. 

The NLO comparison of the complete and NLL expanded results indicates that the latter are 
reliable at x ~ 0.01, but not at x ~ 0.001. The estimated NNLO corrections are large and positive 
around the peak of the distribution, where they amount to as much as 40%. More work is needed to 
arrive at quantitatively reliable NNLO predictions for this and other differential cross sections. It 
is interesting to note, however, that a considerable excess over the NLO results has been observed 



in HERA measurements of charm production including very low values of pj [|23|]. 
5. Summary and Outlook 

We have determined the next-to-next-to-leading logarithmic (NNLL) resummation exponent and 
the one-loop matching function for the dominant gluon coefficient function for the heavy-quark 
structure functions F2 in deep-inelastic lepton-hadron scattering. The results have been used to 
obtain all threshold-enhanced NNLO contributions to this coefficient function, which we have il- 
lustrated for the especially important case of charm production. 

At present, these results provide the only reliable estimate of the NNLO effects at small scales, 
^ lOm^. At larger scales, it may be useful to combine these threshold contributions, the Mellin 
moments (with respect to x) of the large-(§ limits [0] and the leading large-Tj (small-x) logarithms 



[]24[], in order to obtain an all- 77 approximate NNLO coefficient function. 

As an example for less inclusive quantities, we have also presented NNLO threshold estimates 
for the /jy-differential structure function. Also here the accuracy reached in Ref. [ p^ ] needs to be 
improved for quantitatively reliable predictions. Present results indicate considerably larger NNLO 
corrections than for F2 close to the peak of the distribution at rather low values of pj. 
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